The nuclear pseudospin symmetry along an isotopic chain 
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We investigate the isospin dependence of pseudospin symmetry in the chain of tin isotopes (from ^^"Sn until 
^^"Sn). Using a Woods-Saxon parametrization of the nuclear potential for these isotopes we study in detail the 
effect of the vector-isovector p and Coulomb potentials in the energy splittings of neutron and proton pseudospin 
partners in the isotopic chain. We conclude that the realization of nuclear pseudospin symmetry does not change 
considerably with the mass number, and is always favored for neutrons. We also find that the p potential accounts 
for essentially all the pseudospin isospin asymmetry observed and that the Coulomb potential plays a negligible 
role in this asymmetry. This can be explained by the dynamical nature of pseudospin symmetry in nuclei, 
namely the dependence of the pseudospin splittings on the shape of the nuclear mean-field potential. 

PACS numbers: 



I. INTRODUCTION 

Pseudospin was introduced in the late 60's 01|^ to account 
for the quasi-degeneracy of single-nucleon states with quan- 
tum numbers (71, i,j = £+l/2)and(n-l,£+2, j ^ t+312) 
where n, I, and j are the radial, the orbital, and the total an- 
gular momentum quantum numbers, respectively. These lev- 
els have the same "pseudo" orbital angular momentum quan- 
tum number, I — I and "pseudo" spin quantum number, 
s = 1/2. 

The pseudospin symmetry has been analyzed in the frame- 
work of non-relativistic models by several authors |3, 4, 5]. 
The subject was revived in recent years when Ginocchio 
presented a relativistic interpretation for this symmetry show- 
ing that the quantum number ^ is the orbital quantum number 
of the lower component of the Dirac spinor for spherical po- 
tentials. Moreover, he showed that £ is a good quantum num- 
ber in a relativistic theory for the nucleon with scalar S and 
vector V potentials with opposite signs and the same magni- 
tude, i.e., E = 5 + y = 0. Actually, this condition can be 
relaxed to demand that only the derivative of S be zero f?!], 
but, if S goes to zero at infinity, both conditions are equiva- 
lent. 

Unfortunately, neither of these two conditions are met in 
nuclei. Some recent works indicate that the nuclear pseu- 
dospin symmetry has a dynamical character It was shown 
that the almost degeneracy of the pseudopsin doublets not 
only depends in the shape of the nuclear mean fields but 
in fact arises from a significant cancellation among the terms 
that contribute to the energy and not only from the pseudospin 
orbit interaction |9, 10]. The non-perturbative nature of this 
interaction has been discussed in Illnl2il . where it also was 
pointed out the dynamical character of the nuclear pseudospin. 

We have investigated recently in a quantitative way the ef- 
fect of the Coulomb and the vector-isovector p potentials in 
the proton and neutron asymmetry seen in the nuclear pseu- 
dospin fl 3l IT4j] . We analyzed the effect of these potentials 



in each of the terms that contribute to the pseudospin energy 
splitting. We concluded that the isospin asymmetry seen in 
the nuclear pseudospin is also a manifestation of the dynami- 
cal character of this symmetry. 

In this paper we will present new results concerning the 
isospin dependence of the nuclear pseudospin along Sn iso- 
topic chain. We will plot the energy splitting for some pseu- 
dospin doublets as a function of the mass number A, analyz- 
ing separately the neutron and proton levels. To this end, we 
solve the Dirac equation for the nucleons using a parametriza- 
tion for the nuclear potential along this chain which was de- 
veloped by us in a recent work 1 14]. That parametrization was 
obtained by a fit to a sophisticated self-consistent relativistic 
mean field calculation 1 15]. We will show that the pseudospin 
quasi-degeneracy has a small dependence with A in accor- 
dance with 1 16]. 

The symmetry is better realized for neutrons than is for pro- 
tons. Furthermore, as found in previous works [8^ .14] , we 
also show that the main reason for the difference between 
the energy splitting of proton and neutron pseudospin part- 
ners comes from the vector-isovector potential and not from 
the Coulomb potential. 

The paper is organized as follows: in section II we present 
the nuclear potentials as combinations of Woods-Saxon po- 
tentials and the respective parameters as functions of A and 
N — Z. h\ section III we analyze the energy splittings of two 
pseudospin doublets that are close to the Fermi sea, for neu- 
trons and protons, as a function of A from 120 to 170. In 
particular, we look for the effect of the Coulomb and vector- 
isovector potentials in those energy splittings. Finally, we 
draw our conclusions. 

II. DIRAC EQUATION AND PSEUDOSPIN SYMMETRY 

The Dirac equation for a particle of mass m in external 
scalar S and vector V potentials is given by 



[a • p + /3(m + 5) + = e * , 



(1) 



2 



where a and (3 are the usual Dirac matrices. The Hamiltonian 
in Eq. Q is invariant under SU(2) transformations when S = 

or 5* = —V 1 17, 18, 19]. The second case corresponds to 
the realization of pseudospin symmetry. As referred before, 
this symmetry is related to the orbital quantum number of the 
lower component of the Dirac spinor 

If S and V are radial potentials, equation Q gives rise to 
two differential equations for the upper and lower radial wave 
functions. Defining /S. = V — S,Y,^V + S, and the binding 
energy E = e — m, these are 



F' 



1 d 

dr 
1 + K 



Gk, 



A' 



E + 2m- A 



r 

^d_ 

I- K 



(£; + 2m- A)(£:-S)Gk = (2) 



-F, 



dr 

{E + 2m 



A){E- 



E-Y. 
S)F«=0, (3) 



where k is the quantum number defined by 

-(£+1) j=l- 



The term with 1 — k in Eq. (|3} is the pseudospin-orbit term 
lEoll . From that equation one sees that, should it be possible 
to set E' = 0, ^ would be a good quantum number. Since 
the sign of k determines whether one has parallel or antipar- 
allel spin and t — I ~ one sees that pairs of states with 

K = -(^ + 1) and K = ^ + 2 have the same £ = € + 1, the 
quantum numbers of the pseudospin partners. For example, 
for [nsi/2, (n- 1)^3/2] one has [ = 1, for [npg/a, (n-l)/5/2] 
onde has I — 2 etc. Pseudospin symmetry is exact when dou- 
blets with j = I ± s are degenerate. 

III. THE NUCLEAR POTENTIAL FOR THE SN ISOTOPIC 
CHAIN 



The parametrization for Rc and Oc has a natural justification 
in view of the known A^^^ dependence of the nuclear radius. 
The Vop, Rp and ap dependencies on iV — Z are also justified 
since they are proportional to the difference between proton 
and neutron densities. For Vcoui we take the proton electro- 
static potential energy in a uniform spherical charge distribu- 
tion of charge Ze and radius R^ 
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.(12) 



We used in Eq. ([T3 Rc = 1.20^^/^ which has the same A 
dependence of the nuclear radius. 

We show in Fig.^the binding E potential for neutrons and 
protons in the Sn isotopic chain. In the proton case we see 
explicitly the Coulomb barrier produced by the potential il2\ . 
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In a recent work we found a general parametrization in 
terms of Woods-Saxon potentials for the binding potential E 
in the whole Sn isotopic chain. The procedure used to extract 
the parameters is explained in detail in 1 14]. We separated E 
in a central Ec(r), a vector-isovector Vp{r) and a Coulomb 
part Vcoui {r) (only for protons), in the following way 

Eq C I p 



1 + exp[(r 
+Vcoui{r) . 



± 



Rc)/ac\ l + exp[(r 



Rp)/ap] 
(5) 



In Vp the plus sign is for protons and the minus sign for neu- 
trons. The parameters, as functions of A, N and Z, are 



Rc 

Vop 

Ro 



-69.94 MeV 
1.21^1/3 

0.13A^/^ fm 
-[0.12(iV- Z) 
[0.03(7V -Z) + 
[0.007(7V- Z) - 



-f 3.87] MeV 
5.05] fm 
- 0.27] fm . 



(6) 
(7) 
(8) 

(9) 
(10) 
(11) 



-20 



-60 



-80 









/.■'/.'•/,■ 




A 




/# 
/■'■'/'■•■■'' 

'•//■■7' 

,0 


--- 100 
110 
--- 120 

130 

140 

150 

- - 160 
170 






, 1 , 



(b) 



r[fm] 



10 



FIG. 1 : Woods-Saxon potential for neutrons (a) and protons 
(b) along the Sn isotopic chain. 
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IV. ISOSPIN DEPENDENCE OF PSEUDOSPIN 
SYMMETRY 



Using the nuclear potential presented in Eq. (|5Jl we now 
display the energy splitting for the two pseudospin doublets 
[3si/2 — 2^3/2] and [2^5/2 — 1.97/2]- Firstly we present results 
for the neutrons and after for the protons where we need to 
consider also the Coulomb potential. 
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FIG. 2: Energy splittings of the pseudospin pairs [3si/2 — 
2^3/2] and [2o?5/2 — 157/2] as a function of the mass num- 
ber A for neutrons with and without Vp. 

In Fig.|2lis shown the neutron energy splittings for those 
pseudospin partners as A changes. In both cases the energy 
differences increase just slightly as A increases. Since the 
pair [2d5/2 — I37/2] is more bound than [3si/2 — 2^3/2], one 
may notice also that the energy splittings are smaller for levels 
closer to the the Fermi sea. 

In a recent work |8] a systematics was established relating 
the pseudospin splittings to the Woods-Saxon parameters of 
the nuclear potentials. The conclusions were that the energy 
differences become smaller as the magnitude of the well depth 
I E I decreases and as the diffusivity (of both S and A poten- 
tials) increases, but grow when their radius increases. From 
Fig. n we see that in the isotope chain the magnitude of the 
neutron potential well decreases and its diffusivity increases 
as A increases, thus favoring the pseudospin symmetry. How- 
ever, the increase of the radius with A^^^ works in the opposite 
way and thus essentially cancels the previous effects. Thus, 
the dependence with A of pseudospin splittings for neutrons 
is small. 

We now perform the same study for protons. The result 
is presented in Fig.|3] Again, we see that there is almost no 
change of the energy splittings with the mass number A. Also, 
as in the neutron case, the splittings are smaller for higher lev- 
els. For example, for ^4=150 the energy difference for the 
doublet [3si/2, 2^3/2] is 1.19 MeV for protons and is 0.49 
MeV for neutrons, whereas the [2^5 /2 , 157 72] pair has an en- 
ergy difference of 3.58 MeV for protons and of 1.97 MeV for 
neutrons. We see that, for the same value of A, pseudospin 



> 




S E-2d5^2- 1^7/2 

full potential 
@-® without V^^^j 
Q-Q without Vp 

■-■ full potential 
H-s without V^^^ 
without V 




120 130 140 150 160 170 
A 

FIG. 3: Energy splittings of the pseudospin pairs [3s 1/2 — 
2^3/2] and [2(i5/2 — I57/2] as a function of the mass number 
A for protons with and without Vp and Vcoui- 



symmetry is favored for neutrons. 

What is the origin of this experimentally observed isospin 
asymmetry of pseudospin in nuclei? Clearly this must be re- 
lated to both the Vp potential, which is repulsive for neutrons 
and attractive for protons, and the Coulomb potential, a repul- 
sive potential which only exists for protons. 

Since the Coulomb potential only affects protons, we could 
think that it could give an important contribution for this 
asymmetry. However, from Fig. |3l we see that, when the 
Coulomb potential is turned off, the pseudospin splittings re- 
main practically the same. For example, for j4=150, if we re- 
move the Coulomb potential, the difference of energy for the 
doublet [3si/2, 2d3/2] is 1-13 MeV (compare with the value 
of 1.19 MeV given before), whereas for [2^5/2, 157/2] is 3.38 
MeV (3.58 MeV). 

In 1 14] we saw that this fact is due to a cancellation be- 
tween the diverse terms that contribute to the energy of each 
level. In that study, we showed that, while the contribution of 
the Coulomb potential to the pseudospin-orbit term is substan- 
tial, it is cancelled by the contribution of all the other terms. 
Fig-EJshows that, although small, the Coulomb potential con- 
tribution is bigger for the pseudospin splitting of deeper pseu- 
dospin partners. 

In a similar way, we analyze the role of the potential Vp in 
the isospin asymmetry by looking into its effect in the split- 
tings of the neutron and proton pseudospin pairs. In Fig.|2]is 
presented the neutron energy difference of these pairs when 
the potential Vp is excluded. When this is done, the energy 
difference increases, showing that this potential favors pseu- 
dospin symmetry for neutrons, which can be understood by 
the systematics referred to before, since Vp, being positive, 
decreases jSj (note that S < 0), and at the same time makes 
it more diffuse. This can be see in Fig.^for the neutrons: Vp 
is larger as A and N — Z increases and, as a consequence, the 
well potential is less profound and more diffuse. For exam- 
ple, again for ^=150, the energy difference for the neutron 
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doublet [3si/2, 2^3/2] without the p potential is 0.86 MeV 
(0.49 MeV with the full potential), whereas for [2^5/2, 157/2] 
is 2.76 MeV (1.97 MeV). 

In the proton case, shown in Fig.|3]the effect of Vp is the op- 
posite, as could be expected, but is smaller in magnitude than 
in the neutron case. This time becomes bigger when Vp is 
added and less diffuse thus working against the realization of 
the pseudospin symmetry. For v4=150 the energy difference 
for the doublet [851/2,2^3/2] is 0.86 MeV (previously 1.19 
MeV), whereas for [2^5/3, 157/2] is 2.87 MeV (3.58 MeV). 
This analysis allows us to conclude that when the vector- 
isovector potential Vp is excluded the pseudospin asymmetry 
for protons and neutrons almost disappears. Therefore, this 
potential is the main responsible for this asymmetry and the 
Coulomb potential Vcoui does not play a significant role. 

V. CONCLUSIONS 

We have investigated in a quantitative way the isospin de- 
pendence of the nuclear pseudospin along the Sn isotopic 
chain. To do this analysis we performed a mean-field model 
calculation with a parametrization for the nuclear potential in 
this chain developed by us in a recent work 1 14]. We used 
Woods-Saxon potentials to fit the proton and neutron poten- 
tials obtained by a sophisticated self-consistent calculation for 
the Sn isotopes il5H . 

In this general parametrization of the potential we have sep- 
arated explicitly the Vcoui and Vp potentials. These two po- 
tentials are the main source of the isospin dependence of the 



nuclear pseudospin. In order to identify the origin of this de- 
pendence we have analyzed in detail the effect of those poten- 
tials separately in the proton and neutron pseudospin energy 
splittings along the Sn isotopic chain. 

We conclude that the dependence with mass number A of 
the pseudospin symmetry measured by the energy splittings 
is small either for neutrons or for protons along the isotopic 
chain. The effect of the Coulomb barrier is also very small 
and almost negligible for the proton levels close to the Fermi 
sea. The difference seen in nature for the pseudospin energy 
splitting of the neutrons and protons comes essentially from 
the vector-isovector Vp potential. It makes the binding E po- 
tential more diffuse for neutrons than for protons. Thus the en- 
ergy splitting decreases for neutrons and increases for protons 
originating this isospin asymmetry in the nuclear pseudospin. 

Finally, from our analysis we can conclude that, at least 
for tin isotopes, the realization of pseudospin symmetry is al- 
most independent of the neutron content of nuclei in a isotopic 
chain. 
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